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ATTRACT 

This paper considers the problem of designing a !?r»eor controlle '-r sys- 
tem's subject to inequality variance constraints. A quadratic penalty function 
approach is used to yield a linear controller. Both the weights in the quadratic 
penalty function and the locations of sensors and actuators are selected by 
successive approximations to obtain an optimal design which satisfies the input/ 
output variance constraints. The method is applied to NASA’s 64 meter Hoop- 
Column Space Antenna for satellite communications. In addition the solution for 
the control law, the main feature of these results is the systematic determin- 
ation of actuator design requirements which allow the given input/output perfor- 
mance constraints to be satisfied. 


I. INTRODUCTION 

Consider the task of controlling the linear, stochastic syste. 1 

(la) x = Ax + B(u-tw) , xeR 0 , ueR^, yeR* 1 

(lb) y = Cx 


(lc) 


z = Px + v, zeR 


W) 


'x(O) 


E I w(t) 1 = 0, El w(t) I (x T (o), w T (t), v T (t)) = 


,v(t) 


v v (*), 


*0 0 0 

Wfi(t-x) 

0 0 Vfi(t-x). 


W = diag [... ...], V = diag [... 

such that these four control design goals are net: 

(I) E^U) = lir E y. 2 (t) < of, i = 1 k 

t-H» 

Eji 2 (t) = lin E u^ 2 (t) < y.. 2 , i = 1, .... n 
t-H» 
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(II) Only £ < £ sensors are used 


( 2 ) 



= Mx + v 

from the admissible set of £ sensors described from (lc) 
1 


i \ / "i Tx + v i 


(lc) 


z = 



* Mx + v. 


Z £ / \*l X + Y 


(3) 


(4) 


(III) Only m < m actuators are used 

m 

B(u + w) = l b.(u. + w.) 

i=l 1 1 1 

/rom tfee admissible set of n actuators described from (la) 
n 


B(u + w) = l b.(u. + w.) 

i=l 1 1 1 


(IV) The control u(t) is a linear function of the present and 
past measurements z(t), t t. 

Many engineering control design problems can be stated with performance 
constraints of the form (I). For example, large space telescopes are feasible 

2 1/2 

only if the R!!S pointing errors (E v. ) are within certain bounds 
2 1/2 ' 

{Ej/r") ' < a.) so as to achieve diffraction-limited performance (o^) of the 

optics. The designer may also have the freedom to choose from a number of 
different types of sensors and actuators at a number of different locations. The 
locations and the types of actuators (sensors) determine the vectors b. (m. ) in 
(4) and (lc). 1 1 

A straight-forward approach to accommodate the bounded input/output problem 
(I) yields nonlinear controllers [1-2], violating goal (IV). A straight-forward 
approach to accommodate goals (IV) and (I) is to use a penalty function method 
[3-5], minimizing 
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(5) 


t 

t m \ f (llyllq 2+ llu|| R 2 ) dx , ||y|| Q 2 = y T Qy 

*o 


while adjusting Q and R until (I) is satisfied. These successive approximation 
schemes [3-5] presume a fixed measurement/control structure, and hence do not 
satisfy goals (II) and (III). It is Important to unify the treatment of all 
four goals (I-IV) since it has been shown [6-7] they are inherently interdepen- 
dent problems. In particular, for the isolated problems; [6] has shown the 
optimal sensor and actuator selection for LQG problems (5) with fixed (Q,R), and 
[3-5] have adjusted Q and R to satisfy the constrained-variance problem (I) with 
fixed sensors and actuators (i.e. fixed B, H). 

Unfortunately, the optimal answer for the simultaneous solution of both 
problems turns out not to be the juxtaposition of results [6] and [3-5], due to 
the interdependence of the two problems. 

The purpose of this paper is to present a unified treatment of the entire 
problem (I-IV), which we call the Constrained Variance Sensor/Actuator Selection 
(CVSAS) problem. Sectionll describes the approach. Section III gives the formu- 
las for sensor and actuator effectiveness to deal with goals (II)' and (III). 
Section IV presents the numerical algorithm for Iteratively dealing with goal (I). 
Section V gives the algorithm for solving the entire problem (I-IV). Section VI 
illustrates the application to the Hoop-Column Antenna. 


II. APPROACH 

The solution of the problem with inequality constraints (I) is generally 
not unique. To be a bit more specific than statement (I) we define two variations 
of the problem. The first is called the "Constrained-Input Variance" option of 
the CVSAS. In this option the input constraints In (I) are binding and the 
output constraints in (I) are : *lar,ed. 

CIVSAS: The Cons trained -Input Variance. Sensor/Actuator Selection Problem 

Satisfy (II), (III), and with all input-constraints binding > 

-2 2 

(6) u.j E co u 1 * 1, 1*1, ...» m , 

minimize ( recall y^ * c^xj, 

(7a) Vy ' t > 1. 

-2 2 

If however , there is no i for which Vi > 1 then minimize 
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(7b) 



with all input constraints binding (6). 

Definition: The phrase ", minimal achievable output performance" for the CIVSAS 

will man the minimum constraint violation in the sense of the minimum value of 
Vy in (7) with input constraints binding (6). 

The CIVSAS problem is useful when one wishes to determine the best perfor- 
mance achievable for a given power limitation on the input devices (actuators). 
That is, for a given set of y. the CIVSAS finds the minimum achievable output 
performance. 1 

The second variation of the CVSAS problem -is called the Constrained Output 
variance Sensor/ Actuator Selection (COVSAS). ' 

COVSAS: The Constrained-Output Variance, Sensor/ Actuator Selection Problem 
Satisfy (II), (III), and with all output constraints binding 


(8) °rV * * • 

minimize 


i = 1. 


• • • » 


k 


» 


(Sa) 


"u - l V * >• 


(9b) 


If however , there is no 1 for which E^U.. > 1 then minimize 

m p o 

1/ = l E u/ 
u j-j i °° 1 


with all output constraints binding , (8). 

Definition 2: The phrase "minimum achievable input performance " for the COVSAS 

will mean the minimum constraint violation in the sense of (9) , with all output 
constraints binding (8). 

The COVSAS is useful when one wishes to determine the necessary capabilities 
(design requirements) of the actuators In order to achieve the specified output 
performance. That Is, for a given set of o,, the COVSAS finds the minimum 
achievable input performance. 
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III. SENSOR/ ACTUATOR EFFECTIVENESS 

In this Section we temporarily assume that Q and R in (5) are specified 
diagonal matrices Q = diag [... ...], R = diag [... r^ ...], and we wish to 

determine a ranking of the effectiveness of the admissible set of sensors and 
actuators for the LQG problem described by (1) and (5). To help with this task 
a price or "cost" is assigned to each input and output by decomposing the total 
system cost function (5) into contributions from each input and output. This 
task is called "input or output cost analysis" and from [6] we have the results 

m „ k „ m ^ i 
(10) V = I v. u + l v. y = l l/.* + l V* 

i=l 1 1=1 1 1=1 1 i=l 1 

where l/. u , V? t l/. w , is the contribution in V of, respectively, the 1^ 
control u.j, output y., noise w.., or noise v.., and 

(11a) V. u = r-llg-M? i = l, ...,m 

X 

(lib) * qJ |cj | * i = 1»'«..» k 

1 1 1 P-OC 

(11c) V i W = w i i 1 1 b i 1 1 k+L 

(lid) v* = v ii ||f 1 ||f 1 - 1 I 

A 

where P, K, X and L satisfy 

(12a) 0 = PA T + AP - PKVV’P + BWB T , [fj, ...» f £ ] = F = PtfV 1 

(12b) 0 = KA + A T K - KBR -1 B T K + C T QC. l9 v .... gj a 6 T = -KBR" 1 

(12c) 0 = X(A + BG) T + (A + BG)X + FVF T 

(12d) U = L(A - FM) + (A - FM) T L + G T RG 

4»Lk 

The effectiveness of the 1 sensor is measured by 
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(13) 


|/. sens = i/. v 

and the effectiveness of the actuator is measured by 
(14) * V i u - V* 

cpnc 

These terms V. and l/. represent the particular combinations of the input/ 

output costs l/. u , i/- w , l/^ v which are involved in the performance of each sensor 
and actuator. (The distinction here is that the effect of the input can be 
calculated by (/. w , but the effect of an actuator involves both l/^ u and V* since 
the actuator is noisy, and this dependence is accounted for in (14)). To see 

Cane ap+ 

that V. and V. gives the appropriate measure of the effect of deleting 

tn A.U XL X.L. 

the i in sensor or the i w sensor or the l sensor or the i n actuator, refer 
to the numerical work in [7]. 

Two results from [6] add insight into the use of (13), (14). 

Theorem 1, [6,7]: 

For a specified (Q,R)» the optimal value of the LQ6 performance metric (5) 
marmot be reduced by the deletion of any of the admissible sensors z. , 

1 = 1 , .... l. 

Theorem 2, [6,7]: 

For a specified (0,R) the optimal value of the LQG performance metric (5) 
can possibly be reduced by the deletion of some of the admissible actuators 
u.- , i = 1, ...» m. 

cane 

These theorems partially explain why the sensor effectiveness Is a 

much simpler calculation than l/. act . Since the magnitude of the gain on the i** 1 
sensor signal | |f^| | - I l m il Ipp V ii -*• 0 as -*• », an extremely noisy sensor 

simply will not affect the optimal LQG controller. Hence, the effectiveness of 

fU y 

the 1 sensor can be calculated by the Input cost 1/j . Section V will show how 
to use (13) and (14) in the solution of the COVSAS problem. 
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IV. THE COVLQG ALGORITHM 


Now we cite an algorithm (COVLQG) to solve the COVSAS problem under the 
temporary assumption that l = i and m * m. That is, all admissible sensors and 
actuators are used (3 = B and P = M). The COVLQG algorithm will first be stated 
and then its theoretical properties will be discussed. 

The COVLQG algorithm (i.e. the COVSAS with i = m = m) : 

Step A: Compute P from (12a). If o. ^ | | lip > 1 STOP. No solution to 

the COVLQG problem exists. Otherwise initialize 

q^o) = a ." 2 , r^o) = U^ 2 . 


2 

Discussion of Step A: The lower bound on E in an LQG problem is 
E :jf> I l c i! I 2 (from the well known lower bound tr CPC^ on 1/ in (5)), and this 
result is independent of the choice of Q > 0, R > 0. 


Step B: 


Compute 



E 


» U i 


2 



Vi: q.. > 0 


using (11), (12). If a. E V . = 1 V i : q . > 0 and if 

-2 2 11 ' 
y.. ^» u i — ^ V i = 1, ... m, STOP. The COVLQG solution has been 

Discussion of Step B: In the COVLQG option all necessary control effort Is 

2 2 

applied to force the constraints Ejf^ < to be binding. A formal proof 

the stopping criterion of Step B indicates a solution of the COVLQG problem 
given by Theorem 5 of [7]. 


found. 

that 

is 


Step C: Q and R update equations: Let the iteration index be j and set 

q 1 (J+l) = Oj" 2 E oe v. 2 3q i (j) , 1 = 1 k. If (coj 2 )" 1 < q 1 (j+1) < eoT 2 , 

(e < 0 small specified constant) then set q^(j+l) n 0. If 
o i~ 2 E^y.. 2 ■ 1 V 1: q^ > 0, then set r^(j+l) = [y^Eju 2 ]^ 2 ^^), V i: 
y ^” 2 EjiS < 1. For all other i, set r^(j+l) s r^(j). Return to Step B. 

Discussion of Step C: The r.,(j+l) of Step C are clearly adjusted toward the 

' -2 2 

stopping condition of Step B (y^ Ejj^ >_ 1), since a reduction in r^ causes 

Eju* 2 to Increase. The justification for setting q^ = 0 when either q^(j+l) -*■ 0 
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or when q^j+l) -*•<*> is as follows: The tendency of q.. toward zero indicates a 

lack of output controllability due to a degenerate rank of C (rank C < k). In 
this case, the algorithm ceases to attempt the impossible (i.e. to force two 
dependent outputs to arbitrary values) by removing this particular y^ (the least 

critical one as indicated by the smallest q^ -*■ 0) from the cost function by 

setting its coefficient q. ~ 0. Now let rank C = k. The tendency of q. toward 

°° can result only when a stabil izable, detectable system is not output controll- 
able, (even though C = k) and an uncontrollable output converges to a value which 

2 2 

violates its constraint (E jf , , > ). The constraint is violated the smallest 

amount possible since in this case the corresponding q. -*■ <» on successive 

iterations of the update equations. When this condition is determined, such 
y_.'s are removed from the cost function on future iterations (by setting q^ = 0) 

since it now has been established that they cannot be brought within specifica- 

tion.Ej^ 2 < a. 2 . 

A similar algorithm exists for the Constrained Input Variance LQ6 problem 
(CIVLQG) and details are given in [7]. 


V. THE COVSAS ALGORITHM 

The sensor/actuator effectiveness formulas (13), (14) derived in Section III 
and the COVLQG algorithm of Section IV are now integrated to solve the COVSAS 
problem posed in Section'll. 

COVSAS Algorithm : 

- - 2 2 

Step 1. Specify {A,B,C,W,V,*,,m,a ,v >. Run COVLQG algorithm using l 
actuators t m sensors. 

Step 2. Compute y s ® ns f from (13), (14) and rank sensors and 

actuators according to their effectiveness: 


(15a) 


sens 




sens 




sens 


(15b) Vj act > t/ 2 act > ... > V m act 

Delete the sensor and actuator with the lowest effectiveness values 

c AWf 

(/. , I/. , provided such deletion does not cause loss of 
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controllabxlity or observability . * Unless l < £ + 1, reset i xo SL- 

- -? 2 
Inless m < m + 1, reset m to m-1. If o. £ y. =1 V 1 = 1, ...» 

£ 1 1 

and Vi: y. 2 E U . 2 > 1, •£/ [■— 7 y. 2 E u . 2 ]/ • , , \ .. < 

i » i ’ J l £ .£, m i « i J (j+l)iteration 

1 i _ 2 2 

[ T y i E » u i ] (j+l)iteration retw ’ n Ste ? L Otherwise STOP, 
solution to the COVSAS has been found. 


1 . 

k 


A 


Piscussion of Step 2: Nunerica. experience wizh this algorithm suggests that 

more than one sensor and more than one actuator may be deleted on each iteration. 
In fact, tor many cases the sane result can be obtained by reducing t to I and 
m to m on the first iteration. However, this quicker convergence can sometimes 
converge only to suboptimal answers, and the algorithm above is written in its 
most conservative form (deleting only one sensor and/or actuator per iteration) 
where convergence to optimal values is more reliable [7]. 


VI. CONTROL OF A SPACE ANTENNA 

Fig. 1 depicts the Hoop-Column Antenna arrangement for a proposed NASA 
communications satellite. Stationed in a geosynchronous orbit, the objective of 
the antenna control system is to regulate the orientation and focus of the 
satellite antenna relative to its multiple feed horns (at node 10). Table 1 
lists the 24 linear and angular displacements which make up the outputs y. , 

i = 1, .... k, where k = 24. Table 2 lists the 39 admissible sensors and Table 
3 lists the 12 admissible actuators. Note that ARX2 stands for angular rate 
about the x axis at node 2. AX2 stands for angular displacement about axis x at 
node 2. Z10-Z2 stands for a rectilinear displacement between nodes 10 and 2 in 

the z direction. The specifications for the outputs are 0 ^ = 22.8 are seconds 

for i = 1, .... 6, and 0 ^ = .158 mm for i = 7 24. The specifications for 

the inputs u. are y . = 10 dn-cm, i = 1, .... 12. The actuator noise is described 

ii 2 

by W = diag [... ...], = .1 (dy-cm) , V i = 1, ... 12. The sensor noise 

is V = diag [... ...], = 7.615xl0" 7 rad 2 , i = 1,2,3,13,14,15, - 

2.bxl0' 7 m 2 , i = 4, ...» 12, 16, .... 27, V^ * 4.76xl0~ 5 (rad/sec) 2 , i ■ 28, 

..., 39. It is des: ted to limit the number of actuators to 6 » m and the number 
of sensors to 12 = I. The dynamics of the antenna structure were described by 
10 elastic modes and 3 rigid body modes. The square of the frequencies 


^Observability, controllability checks are particularly simple for flexible 
space structures using the tests in [8]. That is, rank tests of matrices 

[B; AB, ... A n_1 B], [C T , A T C T , ... A Tn_1 C T ] can be avoided. 
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2 

id.. , i = 1, .... 10 of the elastic nodes are 

(oij 2 , uj 2 2 , .... a> 10 2 ) = (.40579, 7.2090, 7.2362, 13.2?/, 
44.834, 132.14, 14?. 66, 445.01, 448.69, 775.86) (rad/sec) 2 . 


More complete information for the antenna model nay be found in [7]. 

The results of the COVSAS algorithm applied to the Hoop-Column Antenna are 
summarized in Table 4. The 6 actuators deleted from the admissible set of 
Table 3 are (listed in order of deletion): u g> u g» u iq* u 7* u 4 # ^ 

sensors deleted (in order of deletion) are: z^g, z 3 , Zg, Zj 2 , z, Zj 3 , z 2 , Zj, 

z 24 , z 27 , z 4 , z 5 , z 18 , z 21 , z 3Q , z 3g , z 33 , z 7 , Zg, z 31 , z 23 , k. 20 , z 35 , z 25 , z 22 , 

z^g. Notice that even though the output constraints are still binding the total 

control effort is less using only 6 actuators, (6x5.021 ■ 30.12) than using 12 
actuators (12x3.275 = 39.30 > 30.12). Thus, better performance is possible with 

fewer actuators, since for several actuators the noise effect P. Vi is greater 

than the signal effect V. u in (14) (note thw- negative values of i^. act in Table 

4). 


Perhaps the most important information from t^® C^’^AS is the determine* 
of the minimum achievable actuator specification ' r from Table 5 that all 
of the 24 outputs are held within their design constiainu (a* = 22.8 are secs. 

for angles and a.. = .158 mm for rectilinear displacements) by actuators which 

must be design*, i for the capabilities of TABLE 5. That is, the given output 

specifications, of are possible to meet if is changed (=> actuators are 

redesigned) (from Table 5) to Uj = 73, y 2 = 26, y 3 * 105, y 4 = 26, y g * 32, 

u 6 = 39 * 


VII. CONCLUSIONS 

Presented is an algorithm COVSAS which Integrates the following tasks: 

Selects sensors and actuators from an admissible set. 

Designs a linear feedback controller which satisfies output variance 
constraints. 

Determines actuator design requirements which allow the output variance 
constraints to be satisfied. 
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Numerical properties of the convergence of this algorithm are given for NASA's 
Hoop-Column Antenna. Additional theoretical properties of convergence of this 
algorithm are given in [7], 
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Figure 1: Hoop Column Antenna 
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Table 1: Hoop Column Output Description 

Output # Type Nodal Location Direction 


1 

Inertial Angle 


2 


X 

2 


II 


2 


Y 

3 


II 


2 


Z 

4 

Rel-ative Angle Between 10 

and 2 


X 

5 


II 


M 


Y 

6 


Inertial Angle 

10 


Z 

7 

Relative Linear Disp. 

Between 6 

and 2 


X 

8 


11 


It 


Y 

9 


II 

9 

and 2 


X 

10 


H 


N 


Y 

11 


II 

10 

and 2 


X 

12 


n 


N 


Y 

13 


it 

101 

and 10 


X 

14 


it 


H 


Y 

15 


ii 


II 


Z 

16 


n 

107 

and 10 


X 

17 


ii 


II 


Y 

18 


ii 


II 


Z 

19 


ii 

113 

and 10 


X 

20 


it 


li 


Y 

21 


ii 


H 


Z 

22 


n 

119 

and 10 


X 

23 


it 


II 


V 

i 

24 


it 


II 


z 


Table 2: Hoop- 

Column Sensor Labels 



Sensor 


Sensor 


Sensor 



Number 

Label 

Number 

Label 

Number 

Label 


1 

AX2 

14 

AY 10 

27 

Z119-Z10 

2 

AY2 

15 

AZ10 

28 

ARX2 


3 

AZ2 

16 

X101-X10 

29 

ARY2 


4 

;;6-x2 

17 

Y101-Y10 

30 

ARZ2 


5 

Y6-Y2 

18 

Z101-Z10 

31 

ARX6 


6 

Z6-Z2 

19 

X107-X10 

32 

ARY6 


7 

X9-X2 

20 

Y107-Y10 

33 

ARZ6 


8 

Y9-Y2 

21 

Z107-Z10 

34 

ARX9 


9 

Z9-Z2 

22 

X113-X10 

35 

ARY9 


10 

X10-X2 

23 

Y113-Y10 

36 

ARZ9 


11 

Y10-Y2 

24 

Z113-Z1Q 

37 

ARX10 


12 

Z10-Z2 

25 

X119-X10 

38 

ARY10 


13 

AX10 

26 

Y119-Y10 

39 

ARZ10 
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Table 3: Hoop Column Actuator Description 

Actuator torque abott 

axis at 
Node location 

Uj = T X 2 
u ? = T Y 2 

u 3 s T Z 2 
u 4 = T X 6 
u 5 = TYb 
u 6 = JZ 6 
u ? = T X 9 
u 8 = TY9 
Ug - T Z 9 
u 10= TX10 
u n = T Y 10 
u 12 = T Z 10 
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Table 4: Hoop Column Output Constrained COVSAS Results 


Iteration 

Number 

Identified 
Sensors (l£ en ) 

Identified 

Actuators 

d^ ct ) 

Ave Input 
Value 
(7.6) 

Number of 
Sensors/Actuators 

1 

AZ10(. 00041 16) 
AZ2(. 000397) 
Z6-Z2(0) 
Z9-Z2(0) 
Z10-Z2(0) 

TZ10(-1.362) 
TZ9(— 1.369) 

3.275 

39/12 

2 

AY1(.0C3362) 
AX10(. 003358) 
AY2(. 00226) 

AX2(. 00226) 
Z113-Z10(. 001942) 
Z119-Z10(. 001884) 

TZ6(-2. 1405) 

3.592 

34/10 

3 

X6-X2(. 01457) 
Y6-Y2(. 01455) 
ZlOl-ZlO(.OllO) 
Z1C7-Z10(.0108) 

TX10( -1.2055) 

3.699 

28/9 

4 

ARZ2(. 02844 
ARZ10(. 02232) 
ARZ6(. 02238) 

TX9( -1.2917) 

3.997 

24/8 

5 

X9-X2(.0986) 

Y9-Y2(.0839) 

TX6(-1.4793) 

4.377 

21/7 

6 

ARX6(. 07648) 
ARX2(. 07648) 

— 

C.829 

19/6 

7 

Y107-Y10(. 13395) 
XRY9(.1098) 

— 

4.857 

17/6 

8 

XI 19-X10( . 1557 ) 
XI 13-X10( . 1555) 
X101-X10( .1551) 


4.905 

15/6 

9 

.... 



5.021 

12/6 
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Table 5: Output-constrained Specifications 


E u. 


Output # 

E v. 2 

** l 

Actuator # 

W | 

(minimum achievabl 

1(AX2) 

.015 sec 

1 TX2 

72.91 dn-cm 

2(AY2) 

.015 sec 

2 TY2 

26.145 dn-cm 

3 ( AZ2 ) 

11.588 sec 

3 TZ2 

105.47 dn-cm 

4(AX10-AX2) 

.001 sec 

4 TY6 

26.138 dn-cm 

5(AY10-AY2) 

.001 sec 

5 TY9 

31.750 dn-cm 

6 ( AZ 10) 

12.000 sec 

6 TY10 

38.812 dn-cm 

7(X6-X2) 

.010 ran 



8( Y6-Y2) 

.010 ran 



9(X9-X2) 

.068 ran 



10(Y9-Y2) 

.068 mn 



11(X10-X2) 

.158 ran 



12( Y10-Y2 ) 

.158 ran 



13 ( X101-X10) 

.104 ran 



14( Y101-Y10) 

.158 ran 



15fZ101-Z10) 

.007 ran 



16(X107-X10) 

.158 mn 



17 ( Y107-Y10 ) 

.156 ran 



18( Z107-Z10) 

.008 ran 



19 ( X113-X10) 

. 122 ran 



20( Y113-Y10) 

.158 ran 



21 ( Z1 13-Z10) 

.001 ran 



22(X119-X10) 

.158 ran 



23 ( Y1 19-Y10) 

.091 ran 



24(Z119-Z10) 

.001 ran 




36 



